Exterior - Interior Duality for Discrete Graphs 



00 ' Uzy Smilansky 

' Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 

(SI ' 76100, Israel. 

School of Mathematics, Cardiff University, Cardiff, Wales, UK 



> 

(N 
m 



oo 
O 



E-mail: uzy . smilanskyQweizmann .ac.il 



> 

o 

CO ' Abstract. The Exterior-Interior duality expresses a deep connection between the 

Laplace spectrum in bounded and connected domains in K^, and the scattering 
matrices in the exterior of the domains. Here, this link is extended to the study of the 
spectrum of the discrete Laplacian on finite graphs. For this purpose, two methods 
are devised for associating scattering matrices to the graphs. The Exterior -Interior 
duality is derived for both methods. 



1. Introduction 

The purpose of the present paper is to introduce to the study of the Laplacian on 
O ■ discrete graphs a concept which found several applications in the spectral theory of the 
wave (and the Schrodinger) operator - the exterior-interior duality. It applies for wave 
equations on manifolds which can be partitioned into interior (assumed to be simply 
connected and compact) and exterior domains separated by a boundary. One considers 
now the wave equations in the two domains, subject to the same boundary conditions. 
■ The wave operator restricted to the interior has a pure point spectrum, while in the 
^ ■ exterior the spectrum is continuous and a scattering operator can be defined. The 
restriction of the scattering operator to a particular value of the spectral parameter is 
the unitary (on shell) scattering matrix S{X). Its spectrum is confined to the unit circle. 
The exterior-interior duality asserts that the spectrum of the interior can be identified 
as those values of A for which one of the eigenvalues of the 5* matrix approaches the 
value 1. 

The exterior-interior duality was first discussed in the physics literature in 
[H El [31 H]. It was presented in several versions, which differ in the geometry of the 
space where the scattering operator is defined. Consider e.g., a compact domain i7 C 
and the wave equation subject to Dirichlet boundary conditions. In the exterior, the 
wave equation is subject to the same boundary conditions, and Q is treated as an 
obstacle. The scattering it induces is described in terms of the scattering matrix S{X). 
The exterior-interior duality ensures that the spectrum of the Dirichlet Laplacian in 
the interior is identified by the values of the spectral parameters where an eigenvalue of 
S{X) approaches unity. This was precisely formulated and proved in [5], El [7]. 
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Figure 1. The domain and tlie added auxiliary channels which are used for the 
definition of the scattering matrix. 



The wave-guide version which is closer to the present discussion can be crudely 
described as follows: Consider a connected domain in and impose Dirichlet boundary 
conditions on its boundary dVt which is assumed to be smooth. The Dirichlet Laplacian 
has a discrete spectrum. To define an associated scattering system, choose on the 
boundary two arbitrary points Ai and A2 (see figure ([1])). Draw two parallel lines, 
li through Ai and I2 through A2 so that li and I2 are perpendicular to the chord 
r = (^2,^1). Consider now the part of the boundary dVti which starts at Ai and 
is traversed in the positive mathematical sense as one goes from Ai to A2. The domain 
which is bounded by dVl^ and by the half lines of li and I2 to the left of F is a semi- 
infinite channel, terminated at one of its ends by dVLi. This domain will be denoted by 
L (for "left" ) and its boundaries are drawn by darker lines. Impose Dirichlet conditions 
on the boundary of L and compute the scattering matrix S'l(A). Repeat the same 
construction for dVL^ which is the part of the boundary which completes dVti to the full 
boundary. The scattering matrix Sr{X) in the "right" channel R is defined in the same 
way as above (The domain boundaries are drawn with lighter lines). The spectrum of 
the Laplacian in VL is obtained by solving the secular equation 



This equation expresses the essence of the exterior-interior duality: the spectral 
information in the interior is obtained from the scattering matrix computed for an 
appropriate exterior problem. This version and its variations was studied in the physics 
literature [H El [9] . Only a few examples were studied with mathematical rigor [10] . 
The exterior-interior Duality was discussed and proved for quantum graphs. 



Secular equation of the type (II. ip have several theoretical and practical advantages 
which were used for various purposes in the study of the wave equation. In the sequel, 
analogous secular equations for the spectrum of the discrete Laplacian on graphs will 
be derived, thus extending the interior-exterior duality to the discrete case. There are 



det (/ - S{\)) = ; S{\) = Sl{\)Sr{\) 
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two conceptual obstacles which should be cleared in order to carry out this program: 

1. Given a finite graph with a general connectivity, there is no natural way to convert 
it to a scattering problem analogous to the first version discussed above. 

a. There is no natural definition for a "boundary" on a compact graph, and if a boundary 
is defined, the meaning of boundary conditions has to be explained. 

In the following section the above difficulties will be addressed, and the scattering 
matrices >S'(A) will be constructed in two different settings. The exterior-interior duality 
will be realized by deriving the corresponding secular equations of the form f 1 1.11) for the 
spectrum of the graph. 

A few applications will be presented in section (jlj). 

The following notations and concepts will be used throughout the present work. A 
graph Q is defined in terms of its vertex set V and edge (bond) set £. The cardinality 
of these sets will be denoted byl^=|V|, E = \£\. The graph topology is defined 
by the connectivity (adjacency) matrix C which is labeled by the vertex indices. It 
is a symmetric matrix whose i,j element counts the number of edges connecting the 
vertices i and j. The valency (degree) of a vertex is defined as the number of bonds 
which emanate from it, Vi = Ylj Cij. The diagonal matrix D = diag(f j) appears in the 
definition of the graph Laplacian: 

A = -C + D. (1.2) 

As explained above, and in a complete analogy with the example discussed in the 
introduction, the graph whose spectrum one wishes to compute will be extended so that 
it forms the "interior graph" Q^^^ in a scattering setting. All the quantities which are 
related to the "interior" graph will be denoted by a superscript 0. For simplicity it is 
assumed to be connected and devoid of parallel bonds (Q^j G {0, 1}) or loops {C^'^i = 0). 

2. Scattering on discrete graphs 

In this section scattering on graphs will be defined using two different approaches. The 
first goes along the lines of [13] , the second relies on the graph evolution operator defined 
in [H] . For the sake of completeness, some background material for each of the methods 
will be reviewed before the scattering matrices are presented. 

2.1. Scattering into attached leads 

In this setting, one attaches leads to the graph Q^'^^ and scattering is defined on the 
enlarged graph. A lead graph / is defined as a semi-infinite set of vertices (/, 1), (/, 2), ■ • • 
which are connected linearly. A vertex is identified by a double index (/,«), / denotes 
the lead, and i enumerates the vertex position on the lead. The lead connectivity 
(adjacency) matrix is = w6i^ii6\n-n'\,i , n,n' ^ N+, where w stands for the 

number of parallel bonds which connect successive vertices. (All quantities related to the 
leads will be denoted by the superscript (Lead)). The spectrum of the lead Laplacian 
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^(Lead) corresponding eigenf unctions f = f{i,2), ' ' ')~^ satisfy 

(Aj^''"'^^f = - w(/(/,n+l) + f{l,n-l)) + 2w/(i,„) = Xf(l,n) forn > 1 , 

= -u^/(/,2) + = A/(i,i) forn = l. (2.1) 

The spectrum is continuous and supported on the spectral band A G [0,4w] 
(the conduction band). For any A in the conduction band, there correspond two 
eigenfunctions which can be written as hnear combinations of counter-propagating 



waves: 

For A > 4:W, |^_| > |,^+|. The reason for constructing the leads with w parallel bonds is 
because the conduction band can be made arbitrarily broad. In the present application, 
an appropriate choice of w would be of the order of the mean valency in the interior 
graph, so the spectrum of Q^^'^ falls well within the conduction band. 

A function which satisfies the boundary condition at n = 1 (second line of (12. ip ) is, 

= /(U + ^KA)C , (2.3) 

where 

siW = = ^+ ' I^KA)| = 1 for A G [OAw] . (2.4) 

The lead scattering amplitude s;(A) provides the phase gained by scattering at the end of 
the lead (as long as A is in the conduction band). It plays here the role of the scattering 
matrix which will be defined in the sequel. 

Returning now to the interior graph Q^^^ it is converted into a scattering graph 
by attaching to its vertices semi-infinite leads. At most one lead can be attached to a 
vertex, but not all vertices should be connected to leads. Let C denote the set of leads, 
and L = \C\. The connection of the leads to Q'^^^ is given by the X L "wiring" 
matrix 

I 1 if 7 G V^°^ is connected to / G £ ^, 
= ( otherwise ' <'-^' 

The number of leads which emanate from the vertex i is either or 1, and is denoted 
by di = W^j, (/,!)• Define also the diagonal matrix D = diag((ij) so that 

WW^ = D ; W^W = /(^) , (2.6) 

where I^^^ is the L x L unit matrix. 

The scattering graph Q is the union of Q^^^ and the set of leads C Its vertex set is 
denoted by V and its connectivity matrix C is given by, 

.(0) 



Cl^ if^,jGVW 
Vz,jGV, : C,, = { Ct^^^^^,^^) iileC . (2.7) 

wW,j=(i,i) if i G V(°) and / G £ 
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As is usually done in scattering theory, one attempts to find eigenf unctions f of the 
discrete Laplacian of the scattering graph, subject to the condition that on the leads 
I = 1, ■ ■ ■ , L the wave function consists of counter propagating waves: 

f(i,n) = aiCr' + 6,er ' .n>l. (2.8) 

where a; and hi are the incoming and outgoing amplitudes. They are to be determined 
from the requirement that f is an eigenfunction of the scattering graph Laplacian. It will 
be shown below that this requirement suffices to provide a linear relationship between 
the incoming and outgoing amplitude. The Lx L scattering matrix S^^'^"''^\X) is defined 
as the mapping from the incoming to the outgoing amplitudes: 

b = 5(^""^)(A)a . (2.9) 

To compute ^'•■^'^"'^•'(A), consider the action of the Laplacian on an eigenvector f. 

e : (Af), = - J2 cgV, -wJ2 W^.(U)/(M) + (A + wd,)f, = XU 



V/ G £ : (Af)(;,i) = - w}_^ Wi^i),ifi - ^fm + 2wf^i^,) = A/(,,i) . 

igV(o) 

(Af) - w/(z,„-i) + 2u; = A/(,,„). (2.10) 
The equations for i E V^^^ (first line in (12.101) above), can be put in a concise form: 

(aW + wD- A/^^^"^)) f(^'°') = wWii"^^ . (2.11) 

where /(^'°') is the unit matrix in V^'^'^ dimension, f(^*'''') is the restrictions of f to the 
vertices of the interior graph Q^^^ and f}'^'' is the L dimensional vector with components 
f(i,i), I = 1, ■ ■ ■ , L. For A away from the eigenvalues of A^") +wD the V^^'^ x V^'^^ matrix 
i?(°)(A) is defined as, 

i?(o)(A) = (aW + wD- A/(^*"))"' . (2.12) 

Thus, 

f(v^'°n = ^/?(o)H/fW . (2.13) 

Substituting in the second set of equations in fl2.10p and using (12. Sp . 

{-w^W'^R^^\X)W + {2w - A)/(^)) (a + b) = w;(a e_ + b ^+). (2.14) 

This can be easily brought into the form (12.91) . Using the fact that 2 — X/w = C,- + C,+ 
we get, 

SiLead)^^^ = - {wW^R^^\X)W - e_ /(^))"' {wW^R^''\X)W - e+ . (2.15) 

This is the desired form of the scattering matrix. It has a few important properties. 
i. As long as A is in the conduction band, and are complex conjugate and 
unitary. Since W~^R^'^\X)W is a symmetric real matrix, S^^''°'''-\Xy = S^^'""^\X) and 
g (Lead) (^^^g (Lead) (^^y = j(L) ^ -(.j^g^^ giLead) (^^^ jg Symmetric and unitary matrix. 

a. Once A is outside of the conduction band, the fg^^^ are exponentially increasing 
or decreasing solutions - they are the analogues of the evanescent waves encountered 
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in the study of wave-guides. One of the reason for the introduction of the w parallel 
bonds in the leads was to broaden the conduction band and avoid the spectral domain of 
evanescent waves. However, for the sake of completeness one observes that the scattering 
matrix as defined above can be analytically continued outside of the conduction band 
by using (12. 2p which is valid for any A. The S^^'^"''^'' {X) matrix outside the conduction 
band loses its physical interpretation, and it remains symmetric but is not any more 
unitary. However, it is a well defined object, and can be used in the sequel for any real 
or complex A. In the limit X ^ oo, ^(^^'^^)(A) ^ ^+W^ R^°\X)W ^ (f )^/^^^. 
Hi. At the edges of the conduction band ^±(A = 0) = 1 ; ^±(A = 4w) = —1 . 
Substituting in (12.151) one finds that at the band edges, S^^'^"''^^ = 
iv. The matrix R^^\X) is well defined for A away from the spectrum of A*^"^ + wD. 
Approaching these values does not cause any problem in the definition of S^^'^°''^\X) 
since there S^^^"''^^ = — j(^) . However, for sufficiently large w the sing ularities of /?(°)(A) 
can be separated away from the domain where the spectrum of Q^^^ is supported. 
V. The resonances are defined as the poles of the scattering matrix in the complex A 
plane. They are the solution of the equation 

Zres{X) = det {wW^R'^°\X)W - e_ /(^)) = . (2.16) 

The point A = is not a pole since S'(^'="'^)(A = 0) = -/(^). 

vi. Finally, it might be instructive to note that the matrix R^^\X) is closely related to 
the discrete analogue of the Dirichlet to Neumann map. This can be deduced from the 
following construction: add to each vertex i G V^^-* a new auxiliary vertex i connected 
exclusively to i. (Here we use w = 1 to make the analogy clearer). Write the discrete 
Laplacian for the new graph, and solve (A — AJ)f = 0, where f is a 2V^^^ dimensional 
vector, the first V entries correspond to the original vertices, and the last V entries 
correspond to the auxiliary vertices : f = (f, g)""". Assuming that the values g-^ on the 
auxiliary vertices are given, the entries in f can be expressed as f = R{X)g, where 
R{X) as defined in (12.121) . To emphasize the connection to the Dirichlet to Neumann 
map, define ip = |(g + f) (the "boundary function") and dip = (g — f) (the "normal 
derivative") then, 

= M{X)i) ■ M(A) = 2(/(^"") + i?(A))-^(j(^*°') - R{X)) . (2.17) 
The Dirichlet to Neumann map is defined also in other applications of graph theory, see 

e.g., m- 

2.2. Scattering into dangling bonds (evolution operator approach) 

An alternative construction of a scattering matrix follows naturally from the study of 
the evolution operator [/'^''^(A) of the interior graph. It is a unitary matrix of dimension 
2£^('^) which will be reviewed in the following paragraphs for the sake of completeness. 
(Details can be found in [HI HH]. See also [171 [HI [El EH], and references cited therein). 

Let f = (/i, ■ ■ ■ ; /y(o)) denote an eigenvector of the graph Laplacian A^^^^ (II. 2p . 
corresponding to an eigenvalue A. The bond h connecting the vertices i and j, will be 
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denoted by 6 = {i,j). To each bond b one associates a bond wave function 

ipbi^) = ai, e'^"" + db e-'^"" , x E {±1} (2.18) 
subject to the condition 

M^)=f^ , M-^) = fj- (2.19) 

Consider any vertex indexed by i of degree (valency) Vi, and the bonds {bi,b2, 
which emanate from i. The corresponding bond wave functions have to satisfy three 
requirements in order to form a proper eigenvector of A^^\ 

I. Uniqueness: The value of the eigenvector at the vertex i, ft, computed in terms of 
the bond wave functions is the same for all the bonds emanating from i. The following 
Vi — 1 independent equalities express this requirement. 

ab^ e'^ + db^ e"*^ = ab^ e*^ + db^ e"'^ = . . . = e*^ + e"*^ . (2.20) 

//. f is an eigenvector of A^^^ : At the vertex z, Yl^=i '^^^'^i,jfj ~ ^fi- terms of the 
bond wave functions this reads, 

- J2 k e-*^ + ab, e+^f ] ={\-v,)- J2 [^br. + a,^ e'^f ] . (2.21) 

1=1 ^* m=l 

Together, (12.20p and (12.2ip provide Vi homogeneous linear relations between the 2vi 
coefficients ab^,db^. It is convenient to introduce at this point the following notation: 
Let d = {j,i), i,j G V'-'^-* denote a directed bond pointing from i to j. Then o{d) {t{d)) 
stand for its origin i (terminus j) vertices, respectively. From now on, the amplitudes 
Ofe^ and db^ which refer to propagating on the bond bm in opposite directions will be 
denoted by ad and where d is the directed bond inverse to d. Using f l2.20p and f l2.2ip . 
the outgoing coefficients are expressed in terms of the incoming ones, 

d' : t{d')=i 



where, 



cr 



(0 



I'W = ^ ^Id' - V. A. • (2.23) 



Vil- id 



The vertex scattering matrices a^^\X) are the main building blocks of the present 
approach. A straight- forward computation shows that for real A the vertex scattering 
matrices are unitary and symmetric matrices. 

///. Consistency : The linear relation between the incoming and the outgoing coefficients 
f l2.22p must be satisfied simultaneously at all the vertices. However, a directed bond 
{i,j) when observed from the vertex j is outgoing, while when observed from i it is 
incoming. This consistency requirement is implemented by introducing the Evolution 
Operator Ud'.dW in the 2E^^^ dimensional space of directed bonds, 

Ud',dW = kd),oid')4fW- (2.24) 
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The evolution operator is unitary U = jC^^^*"-') for A G M due to the unitarity of 
its constituents a^^\ Denoting by a the 2£'^°) dimensional vector of the directed bonds 
coefficients ad defined above, the consistency requirement reduces to, 

f/(A)a = a. (2.25) 

This can be satisfied only for those values of A for which 

detfj(2^''")-[/(A)) = 0. (2.26) 



This result can be interpreted in the following way. The evolution operator f/(A) is 
defined for any A, and it maps the 2E^^'^ dimensional vector space of amplitudes a to 
itself. The spectrum of A'^'^) is identified as those values of A for which there exist vectors 
which are stationary under the action of the mapping. 

To construct a scattering operator for which Q^^^ is the "interior", add to the vertex 
set V^^-* another set of vertices denoted by C with L = \C\. Connect £ to a subset of 
V*-°^ such that each / G £ is connected to a single vertex in V*-°\ while a vertex in V^^^ 
can be connected to several vertices in C The vertices in V*-'''' connected to C will be 
referred to as boundary vertices. The V^^^ x L "wiring" matrix is defined similarly to 
the previous definitions (see (12.51) ). 

J 1 if J G V^°^ is connected to / G £ , , 

Wj I = \ , ■ (2.27) 

otherwise 

The number of "dangling bonds" which emanate from the vertex iis di = Yliec 
di can take any integer value or 0. The diagonal matrix with elements di will be denoted 
by D and the identities (12. 6p hold. The new graph Q is of cardinality V = V^^^ + L and 
it consists of the interior graph with L dangling bonds attached. 
The evolution operator for Q can be written in block form as 

/ s ^] \ 

U= p(°"*) , (2.28) 

^ 1]*^ p(*") y 

which is a 2(E^^''+L)x2{E^^^+L) unitary matrix arranged in the following way. The first 
rows and columns are labeled by the indices of the directed bonds which belong to 
the interior graph The last 2L rows and columns are labeled by the indices of the 
directed dangling bond. The first L correspond to incoming bonds pointing from £ to 
V^"-*, the other L indices correspond to outgoing bonds from V*-''-* to C. The 2E^'^^ x 2£'(°) 
upper left block, denoted by S is obtained by modifying the evolution operator for 
the interior graph: the vertex scattering matrices ( 12. 23^ are modified by replacing the 
original valency Vi by the modified valencies Vi = Vi+di. The rectangular 2E^^^xL matrix 
Q consists of elements of the vertex scattering matrices between incoming dangling bonds 
and directed bonds in Q^^^ which emanate from boundary vertices in V'^''''. ^l^"^ provides 
the scattering matrix elements for the time reversed transitions, from the directed bonds 
in Q^^^ to the outgoing dangling bonds. The L x L matrix p(°"*) consists of matrix 
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elements from outgoing to incoming dangling bonds which occur at the dangling vertices 
C Using fl2.23p with vi = 1, we get 



(out) 



l + i{l-\) 



So{d'),t{d)Wo{d'),t{d) 



(2.29) 



1-2(1 -A) 

The L X L matrix p*^*"-* consists of matrix elements from incoming to outgoing dangling 
bonds which occur at the boundary vertices. They are computed using (12.231) with 
Vi = Vi + di. Note: the number of dangling bonds connected to the same vertex can 
exceed 1, and in these cases scattering to other dangling bonds may occur, as long 
as both the incoming and outgoing bonds are connected to the same vertex of Q^^^ . 
The four zero blocks correspond to transitions between bonds which do not follow each 
other. In particular, the second row of blocks in (12.281) has two zero entries, expressing 
the fact that ingoing dangling bonds cannot follow neither interior bonds nor any ingoing 
dangling bonds. 

The unitarity of If ensures the following relations between its components: 

/(^) (2.30) 




a. 
b. 



c. 
d. 



The evolution operator U acts in the 2E^^^ + 2L dimensional space of amplitude 
vectors. Denote the 2E^^^ amplitudes which refer to directed bonds in Q^^^ by a 
and the amplitudes associated with the incoming/outgoing dangling bonds by b*^^/+) 
respectively. The consistency condition (12.251) which determines the spectrum of Q reads 
now, 



(2.31) 



To define a scattering matrix, one has to obliterate the requirement that scattering 
occurs at the dangling vertices. This is achieved by replacing p(°"*) by a block in 
U. By doing so, the rank of U is reduced by L. Thus, the resulting equations do not 
determine a spectrum. Rather, they can be solved for every value of A, and yield a 
linear relation between the amplitudes of the outgoing and incoming dangling bonds. 





Q 









\ 


/a 


\ 








p{out) 






) 


= b(-) 






p{in) 





) 




; 


\ b(+) 


/ 



b(+) 

5(^)(A) 



(^pM + (j(2i?(°)) _ s)-i^]) b(-) = S(^)(A)b 

p(^")(A) + n'^iX) (J^^^'"') - S(A))-ifi(A) . 



(2.32) 



The L X L matrix 5''''^''(A) (The superscript "D" stands for Dangling) provides 
the scattering amplitudes between incoming and outgoing dangling bonds. It is 
meromorphic in A, has a finite number of poles in the lower half plane, and is unitary 
on the real A axis. These statements can be easily checked by using the explicit form of 
the matrix elements of f/(A) and the identities (I2.30p . The main advantage of 5''''^-'(A) 
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over the matrix 5'^^^'*°') (A) defined in fl2.15p is that it depends solely on properties of the 
interior graph, and its definition is free from restrictions or special features due to the 
properties of the leads. 

3. The exterior-interior duality 

In the present section, two versions of the exterior interior duality will be formulated 
corresponding to the two ways by which scattering was defined above. The conversion 
of a graph into a scattering graph can be done in many ways - the vertices connected 
to leads or to dangling bonds can be chosen arbitrarily. It will be shown here, that 
any scattering matrix can be used to extract information about the spectrum of the 
"interior" graph. However, there exists a unique construction which provides the entire 
spectrum of the interior. Namely, when all the vertices are connected to leads (in the 
first setting) or dangling bonds (one per vertex, in the second setting): L = V^^'^ and 
Wi^i = /(^'°'). Under these conditions the secular equations 

det(j(^*'") - S{X)) = 0, with S = OT S = 5^^) (3.1) 

provide the complete spectrum of the interior graph with multiplicities. 

3.1. The exterior-interior duality for scattering to leads 

Equations (12.101) form the basis of this discussion. They can be solved for any value of A 
and they provide the building blocks for constructing the scattering matrix 5'^'''^'^'^'' (A). 
Consider the subset ( 12. lip of these equations. If for a given A the corresponding vectors 
f(^*°') and f^-^) which solve (I2.10p also satisfy 

then, f(^''°^) is an eigenvector of the interior graph Laplacian with an eigenvalue A. If 
the leads are not connected to all the vertices, that is L < V^'^\ the kernel of D is not 
empty. Hence, it is possible that the condition L)f(^*'''') = l^^f^^^-* is satisfied trivially by 
f(^*''') which is in the right kernel of D, and f^^^^ = 0. Such cases have to be excluded 
from the following derivation and therefore, in general, the equation to be derived below 
provides only a sufficient spectral condition. However, if one connects every vertex to a 
lead, this problem does not arise, and a proper secular equation can be derived, as will 
be shown in the sequel. 

As long as A is away from the sing ularities of R^°\X), one can proceed and obtain 
from Df^^'"') = lyff ) 

(wDi?(°) - /(^'°^))iyfj^) = . (3.3) 

Using (12. 6p one finds that (13. 3p is equivalent to requiring a non trivial solution for the 
equation 

W {wW^R^''\X)W - /(^)) ff ^ = . (3.4) 
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Thus, it is sufficient but not necessary for A to be in the spectrum if it is a zero of the 
function 

zr{X) = det {wW^ R^°\X)W - J(^)) . (3.5) 
However, when all the vertices are connected to leads, W = /(^) = /(^*°') , D = /(^*°') 
and i?(°)(A) = (^A(°) - (A - m;)/(^*°')) . Substituting in ([33]), one gets 

, detfAW-AJ(^'"^)) 

--(^^^(-^^^ det(A(o)-(A-.)/(A) - ^'-'^ 
At this point the importance of the parameter w becomes clear: The spectrum of the 
Laplacian coincides with the zeros of -2_r(A) only if there are no accidental overlaps with 
thepolesof i?W(A). Thus, one should choose w so that the spectrum of (A*^") + tt;/(^'°^)) 
exceeds the maximum of the spectrum of A'-^-'. For v regular graphs taking w > 2v is 
sufficient. Using (I2.15P one writes 

wW^R^^\X)W = (/(^) + S^'"""'\X)y^ (e+ /^^^ + 5(^^"'')(A)). (3.7) 

Substituting in (13. 5p yields the secular equation 

^^^(A) = det (J(^) - sr'(A)5(^^"'')(A)) = det (J^^^ - ^-5^^'"''^^)) = . (3.8) 

This is the desired result, since it determines the spectrum of the interior graph in terms 
of the scattering matrix S^^'^'^'^\X) and the free lead scattering matrix si{X) (12.41) . It is 
analogous both in form and in content to the exterior-interior secular equation which 
was introduced in the introduction (ll.ip . 

3.2. Exterior -Interior duality for scattering to dangling bonds 

The scattering matrix 5'''^^(A) will be used to compute the spectrum of the interior 
graph For this purpose, consider the action of the Laplacian of the extended graph 
Q. Denote by fj the components of an eigenvector on a vertex j G V^^\ and by gk the 
component of the same eigenvector on a dangling vertex k & C Then, 

L 

V z : - 5^ a,,/, - ^ W.^igi + {v, + di)fi = Xfi . (3.9) 

A is an eigenvalue of the interior graph, corresponding to an eigenvector f = 
(/i,---,/y(o)), if for every ieV(°) 

L 

J2Wi,i9i = d^fi , (3.10) 
1=1 

and f 7^ 0. Consider a single boundary vertex j with dj ^ 0. Expressing the 
eigenfunction on the dangling bonds in terms of the incoming and outgoing amplitudes 

Wk : W,,k = 1 , fj = bi-^ e-^f + bi+^ e+*f , 

9k = bi'^ e+^f + e-^f . (3.11) 
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The upper line above imposes dj — 1 linear relations between the amplitudes. 
Substituting (13.111) in (13.101) results in one more linear relation. Thus, one gets for 
each j exactly dj linear relations between the incoming and outgoing amplitudes on the 
dangling bonds connected to the interior vertex j. Solving them for each j separately 
and combining the individual relationships gives: 

bW = 5(0)b(-) . ^(°)= 5^ . (3.12) 

This A independent matrix of dimension L x L is bloc diagonal, with blocs of dimensions 
dj. It is easy to check that each block (and hence the entire matrix) is unitary and 
symmetric. When only one dangling bond is connected to a vertex, dj = 1 and the 
corresponding Sj^l = 1. 

The requirement (13.120 has to be combined now with the relation (12.320 which is 
valid for all A. This implies that A is an eigenvalue of the Laplacian on the interior 
graph if there exists a non trivial solution for the equation 

b(+) = 5(^)(A)(5(°))^bW , (3.13) 

in other words, if A is a zero of the secular function: 

zd{X) = det (j(^) - 5(^)(A) . (3.14) 

The derivation above follows from the requirement that (I3.10p is satisfied for f 7^ 0. 
This can be guaranteed only when a single dangling bond is connected to each vertex. 
In this case, 

^(0) ^ j(Vo) and the secular equation (13111 expresses the exterior - interior 
duality in the present setting. 

4. Examples and Applications 

In this section a few examples and applications which will illustrate the exterior - interior 
duality will be discussed. 

4-1- Scattering on a single lead connected to a graph 

Consider a single lead connected to one of the vertices of Q^'^^ which is denoted by i = 1. 
Following (I2:T21[2J[5D . 

G'i,i(A) 



Ri^i{X) = W'R{X)W 



l + wGi,i(A) 



G,,(A) = ((AW - A/(''))-),^, = Y: fry- (4-1) 

r=l ^ 

Clearly G{X) is the Green function (resolvent) of the graph Laplacian. The one 
dimensional scattering matrix reads, 

SiLead)^^^^_ ^Rl,l-^+ , (4.2) 
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and the secular equation: 



1 + 1 



(4.3) 



w;i?i,i(A) -e- l+u»Gi,i(A) 



This function has no poles for A in the conduction band, and it vanishes at the spectrum 
of A^^^) (= the poles of the Green function), provided that the corresponding residues 
do not vanish. This cannot be guaranteed since the component of an eigenvector can 
vanish on any number of vertices. Such eventualities cannot be excluded and therefore 
the proper secular equation has to be constructed by connecting all the vertices to leads. 

4-2. Gomposition 

Consider two graphs Q^^^ and Q^'^^ with vertex sets V*-*^ of cardinality V^^^ and 
connectivity matrices C^*^ ,2 = 1,2. Connect the two graphs by an arbitrary number of 
leads (j, i), i E V^^\ j E V*-^-* but avoid parallel edges. Denote the number of connecting 
edges by L. The connection between the two graphs is given by the matrix: 



Consider the graph Q^^^ say, and regard for the moments the bonds which connect it to 
as dangling bonds. The corresponding wiring matrix is C, and one can write an 
L dimensional scattering matrix 5'^^-'(A) for this system, by using the "dangling bond" 
construction. Similarly, the scattering matrix S^'^\\) which is also L dimensional can 
be written for Q'^'^\ Denote by a and b the L dimensional vectors of amplitudes for 
incoming waves to V*-^-*. The consistency requirement of section (12. 2p requires that 



so that a sufficient condition for A to be in the spectrum of Q^^^ is that it is a solution 
of the equation 



This is the analogue of the example given in the introduction of the exterior - interior 
duality in M^. 

As a simple corollary of this result, one can study the effect of adding a single vertex 
m of valency Vm to a given graph of cardinality V. Assuming that the new vertex is 
connected to the vertices 1, - ■ ■ ,Vm, one computes the corresponding Vm x Vm scattering 
matrix S^^'^ from the original graph. S^^^ plays the role of S^^\\) above. For S^'^\X) 
one uses the vertex scattering matrix a^"^^ f l2.23p . The resulting secular equation reads 




(4.5) 



b = S«(A)a ; a = S(2)(A)b ^ h = S^^\X)S^^\X)h 



(4.6) 



det(j(^)-^«(A)5^'^(A)) = . 



(4.7) 



now 



(4.8) 
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The same problem can be addressed in a different way. With the same construction 
in mind, one can compute the -R(A) matrix (with w = 1). Connecting together all the 
Vm dangling bonds to a single vertex, and setting the value of the wave function to fm 
at the vertex m, one finds for the values of the function on the m connected vertices: 

/, = /^J]i?(A),,fe. (4.9) 

k=l 

At the same time, to be an eigenvalue f must satisfy 

Vm 

Thus, the secular equation for the spectrum of the enlarged graph reads 

Vm 

k,j=l 

It is hoped that the few simple illustrations given above illustrate the advantages and 
potential applications of the exterior-interior duality in the present context. 
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